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ABSTRACT 

Two multifrequency campaigns were carried out on OJ287 in 2005: in April when it 
was in its pre-outburst state, and in November, during the main 12 yr cycle outburst. 
The wavelength coverage was from radio to X-rays. In the optical-to-UV range the 
differential spectrum between the observations has a bremsstrahlung spectral shape, 
consistent with gas at 3 x 10 5 K temperature. Our result supports the hydrogen column 
density of the OJ287 host galaxy of ~ 9.3 x 10 20 cm -2 , the average value found by 
Gosh & Soundararajaperumal. The 3 x 1Q 5 K bremsstrahlung radiation was predicted 
in the binary black hole model of OJ287, and it arises from a hot bubble of gas 
which is torn off the accretion disc by the impact of the secondary. As this radiation 
is not Doppler boosted, the brightness of the outburst provides an estimate for the 
mass of the secondary black hole, ~ 1.4 x 10 8 solar mass. In order to estimate the 
mass of the primary black hole, we ask what is the minimum mass ratio in a binary 
system which allows the stability of the accretion disc. By using particle simulations, 
we find that the ratio is ~ 1.3 x 10 2 . This makes the minimum mass of the primary 
~ 1.8 x 10 10 solar mass, in agreement with the mass determined from the orbit solution, 
1.84 x 10 10 solar mass. With this mass value and the measured K- magnitude of the 
bulge of the host galaxy of OJ287, the system lies almost exactly on the previously 
established correlation in the black hole mass vs. K-magnitude diagramme. It supports 
the extension of this correlation to brighter magnitudes and to more massive black 
holes than has been done previously. 
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1 INTRODUCTION 

There is a fundamental relation between the central black 
hole mass and bulge luminosity in galaxies (Magorrian et al. 
1998, Ferrarese & Merritt 2000, Giiltekin et al. 2009, Kor- 
mendy and Bender 2011a,b). In the study of the tightness 
of this correlation one should use the most accurate mass 
values based on the dynamics of bodies orbiting the central 
black hole. So far there are not many such mass determi- 
nations. The mass of the black hole in the Galactic centre 
gives a fixed point at the low end of the scale (Ghez et al. 
2003, Genzel et al. 2010), the masses of megamaser sources 
fall in the middle of the range (e.g. NGC 4258, Miyoshi et al. 
1995, Kuo et al. 2011), while the mass of the OJ287 primary 
gives the only point at the high mass end (Valtonen et al. 
2010a). In addition, there are plenty of data based on the 
velocity dispersions of galactic centres (e.g M87, Gebhardt 
& Thomas 2009), reverberation mapping measurements and 
finally much recent data based on emission line widths in 
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quasars (for a review, see Ferrarese & Ford 2005). These 
may be viewed as giving reliable but less accurate mass val- 
ues (Valluri et al. 2004, Denney et al. 2009, 2010, Kelly et 
al. 2010). Even though the single data point at the high 
mass end is very accurate, with one percent uncertainty, one 
would like to have independent evidence that it is at least 
approximately correct. Then one could use the dynamical 
mass with greater confidence in correlation diagrams. Thus 
the primary objective in this paper is to confirm the value 
of the central black hole mass in OJ287. 

OJ287 has a historical optical light curve extending 
over 120 yr from late 1800's up to date. It shows a promi- 
nent twelve year cycle (Sillanpaa et al. 1988, Valtonen et al. 
2006a), with an additional two peak structure with 1 - 2 yr 
peak separation (Valtonen et al. 2010b). It was suggested by 
Sillanpaa et al. (1988) that the 12 yr structure arises from 
a black hole binary system where a secondary perturbs the 
accretion disc of the primary at regular intervals and causes 
indirectly increased emission in the jet of OJ287. The jet is 
lined up with the observer at a small viewing angle which 
further enforces the amplitude of brightness variations via 
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Doppler boosting (Lahteenmaki & Valtaoja 1999, Jorstad et 
al. 2005). It is extremely unlikely (x 2 probability less than 
10 -8 ) that that a single black hole model would produce 
such good light curve predictions as the binary model has 
done in the past 14 years (Valtonen et al. 2011). 

This model lead to the prediction that the next major 
outburst after 1983 should take place in the autumn of 1994. 
The outburst came as scheduled (Sillanpaa, et al. 1996). At 
this point several more detailed binary black hole models 
were put forward, with differing predictions for the next ma- 
jor outburst (see Kidger 2000). These were in two general 
categories: exactly periodic (non-precessing) models which 
predicted the main outburst in the autumn of 2006, and pro- 
cessing models which gave the outburst time a year earlier 
(Lehto & Valtonen 1996, Sundelius et al. 1997). There are 
two reasons for this. First, in the latter models the outbursts 
are related to the impacts of the secondary on the accretion 
disc of the primary. With forward precession, as in General 
Relativity, the impact times happen earlier than in a non- 
precessing periodic model. This makes about ~ 6 month 
contribution to the advance of the outburst time (Lehto & 
Valtonen 1996). The second contribution comes from the 
fact that the accretion disc is lifted up by the approaching 
secondary (Sundelius et al. 1996, Ivanov et al. 1998) and the 
impact happens earlier than expected in a rigid disc model 
by another ~ 6 months. 

In order to vertify the correctness, or otherwise, of 
these two clearly different models, a major monitoring cam- 
paign was set up, with frequent observations starting in 
early 2005. Not only was the timing important, but also 
the nature of the radiation at the outburst peak. Most pe- 
riodic models were predicting an increase of radiation from 
the jet, i.e. the outburst should have a power-law like syn- 
chrotron spectrum. In contrast, the precessing binary model 
predicted that the radiation is bremsstrahlung from hot gas 
(1 — 4 x lO 5 ^, the exact value is calculated below) which has 
been torn off the accretion disc. Thus, it was important to 
study the spectral energy distribution during the outburst, 
and to compare it with the spectral distribution outside the 
outburst time. We had two periods of XMM-Newton obser- 
vations, one in April 2005 prior to the outburst, and another 
one in November 2005 when the outburst was expected in 
the precessing model. Preliminary reports of these observa- 
tions have been given in Ciprini et al. (2007), Ciprini and 
Rizzi (2009) and Ciprini (2010). 

As to the timing, the outburst followed the prediction 
of the precessing model (Valtonen et al. 2006b, Valtonen 
et al. 2008). The primary outburst came in 2005 Octo- 
ber/November while OJ287 was in an exceptionally quiet 
state in the autumn of 2006 when the periodic peak bright- 
ness should have happened. Thus the result was unusually 
clear. However, it is also crucial to verify the second hy- 
potheses of the precessing binary model, i.e. to see if the 
spectrum agrees with bremsstrahlung radiation. A major 
problem here would be the contribution from the "normal" 
background activity which produces secondary peaks of radi- 
ation. Luckily the second XMM-Newton period fell between 
two secondary peaks, and thus our measurements should tell 
about the nature of the primary outburst (Figure 1) . 

In this paper we first discuss the differential energy dis- 
tribution between the two XMM-Newton observing periods, 
and confirm that the excess has likely bremsstrahlung na- 
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Figure 1. The optical light curve of OJ287 during the 2005 out- 
burst. Data points are based on Valtonen et al. (2008, 2009). 



ture. Since the radiation from the disc impact is not Doppler 
boosted, the strength of the outburst is directly linked to 
the mass of the secondary. We then ask what is the mini- 
mum mass of the primary which provides the stability for 
the accretion disc. This gives us the minimum value for the 
primary mass. Finally we compare this with the dynamical 
mass determination, and after being satisfied that the two 
methods give consistent results, produce a slightly modified 
black hole mass - K-magnitude correlation diagramme with 
OJ287 added to the older data. 



2 OUTBURST SPECTRAL ENERGY 
DISTRIBUTION 

The two X-ray observations of OJ287 by XMM-Newton in 
2005, April 12 and November 3-4 were performed during a 
ground-based coordinated multiwavelength (MW) campaign 
including both simultaneous radio/mm, near-IR and optical 
observations (provided mainly by the WEBT consortiurrQ) 
and longer term data, obtained within monitoring programs 
by the ENIGMA network facilities, and other participating 
observatories (see, e.g. Ciprini et al. 2007, Ciprini 2010). 

In Figure 2 we plot the observations of the 2005 cam- 
paigns. The observations relating to the two instances of 
time are colour-coded and use different symbols. In addi- 
tion, a collection of data obtained at other times are shown 
by grey open symbols. 

For the 2005 April 12 simultaneous spectral energy dis- 
tribution (SED), the spectrum from radio to X-rays followed 
nicely a synchrotron self-Compton (SSC) model with suit- 
ably chosen parameter values (see the discussion in Ciprini 
et al. 20010). In contrast, the November 3-4 spectrum was 
quite different. The flux came up prominently in the op- 
tical/UV region, but in the radio or hard X-rays the flux 
remained at the pre-outburst level, making the single zone 
SSC model unlikely. 

In Figure 3 we look at the optical/UV region in more 



1 www.oato.inaf.it/blazars/webt/ 

2 www.lsw.uni-heidelberg.de/projects/enigma/ 

3 the names of the 91 WEBT+ENIGMA participants in the mul- 
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Figure 2. The SED of OJ287 during the 2005 outburst in Novem- 
ber 3-4 (blue/dark filled triangle and square symbols), and the 
spectrum of the pre-burst state in 2005 April 12 (red/lighter filled 
diamond and circles symbols) with a pure one-zone SSC fit at- 
tempt (red/light line)). Other archival and not simultaneous data 
(light gray and open pentagon symbols and lines) are reported in 
the plot. 



detail, and plot the April 2005 spectrum as well as the dif- 
ference spectrum (November - April) in Figure 4. The lower 
points (filled squares) are corrected for the Galactic absorp- 
tion (Schlegel et al. 1998, Draine 2003), the upper points 
(crossed squares) are also corrected for the absorption in 
the OJ 287 system. The latter correction is important as 
the Wl channel, represented by the last point on the right 
hand side, lies right on the 2175 A absorption feature of the 
Galactic absortion law (Fitzpatrick 1999). Thus this point 
is a sensitive probe of the interstellar medium of the host 
galaxy of OJ 287. Other points are also affected, but to a 
lesser degree. The hydrogen column density of OJ 287 that 
brings the Wl channel point to the smooth continuation of 
the data from other bands is ~ 9.3 x 10 20 cm -2 . Together 
with the hydrogen column density of 3 x 10 20 cm~ 2 for our 
Galaxy, they add up to the average of the two measured val- 
ues by Gosh & Soundararajaperumal (1995). However, we 
should note that this result is based on the Galactic red- 
dening coefficient Rv = 3.1. If this coefficient is smaller, 
then also the hydrogen column density needs to be smaller 
in order arrive at a fit similar to the one shown in Figure 4. 
For example, Ay = 1.1 implies that the hydrogen column 
density is only half of the above mentioned value. 

The April data shown in Figures 3 and 4 agree well with 
previous measurements from far infrared to UV . The spec- 
trum gradually steepens from the optically thin spectrum 
F v ~ v 
F v ~ i 



• in far infrared (Impey & Neugebauer 1988) to 
in optical (Holmes et al. 1984, Hagen-Thorn 
et al. 1998, D'Arcangelo et al. 2009) to F u ~ v' 17 in ul- 
traviolet (Worrall et al. 1982, Maraschi et al. 1983, Gosh & 
Soundararajaperumal 1995, Pursimo et al. 2000). The steep- 
ening is exactly in tune with the increasing galactic extinc- 
tion from far infrared to ultraviolet in the OJ287 host galaxy, 
and agrees in magnitude with the extinction by the hydro- 
gen column density around 9.3 x 10 20 cm -2 . The standard 



tiwavelength campaign are found in the author list of the Ciprini 
et al. (2007) paper 
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Figure 3. Near-IR, optical, near-UV simultaneous observations 
of OJ287 during the 2005 outburst in November 3-4 (blue/dark, 
top triangle symbols dashed line), compared with the observa- 
tions in its pre-outburst state in 2005 April 12 (red/light bottom 
diamond symbols and dashed line) . Data points are based on the 
same campaign thanks to XMM-Newton UVOT telescope obser- 
vations and other ground-based (WEBT and ENIGMA observing 
campaigns) observatories, from K-band to W2-band. 



extinction law (see e.g. Draine 2003, Fig. 9) for this hydrogen 
column density and shifted to the redshift of OJ 287 z = 0.3 
steepens the spectrum in near infrared by 0.15 units, in op- 
tical by 0.45 units and in UV between A1100 A and A3900A 
by 0.9 units. Other causes for this spectral steepening are 
also possible (Rieke & Kinman 1974), but the coincidence 
with the galactic extinction law is remarkable. 

We have considered two fits to the corrected points: a 
power-law fit with spectral index of a = —0.4 ± 0.17 and a 
bremsstrahlung fit with temperature T = 3^® x lO 5 ^. The 
power-law fit would imply a spectrum which is flatter than 
the normal spectrum by Aa ~ 0.4 units. Such flattening 
of the spectrum at high brightness has been previously ob- 
served in other outbursts (Zheng et al. 2008, Tavecchio et al. 
2010) even though by a smaller amount. In normal flares of 
OJ287 the whole spectrum up to optical-UV is raised, while 
the inverse Compton scattering in the relativistic jet raises 
the X-ray flux correspondingly (Seta et al. 2009). No evi- 
dence of the inverse Compton contribution to the difference 
spectrum is seen in the case of the 2005 November outburst 
(Figure 2). 

The total luminosity of OJ287 from radio to UV would 
be about 2 x 10 46 erg/s if its radiation is isotropic (Worrall et 
al. 1982). However, the radiation is strongly beamed to our 
direction (Sitko and Junkkarinen 1985, Jorstad et al. 2005), 
and thus the total luminosity, integrated over all directions, 
may be only of the order 2 x 10 45 erg/s. It implies a rather 
low accretion rate of ~ 10 -3 of the Eddington rate M-Edd 
for the primary mass of M = 1.84 x 10 10 solar mass. 

In the binary model the secondary impacts the accretion 
disc at the distance of 0.065 pc from the primary in 2005 
(Valtonen 2007). Ivanov et al. (1998) estimate that the peak 
bremsstrahlung luminosity arising from the impact is about 
equal to the Eddington luminosity of the impacting black 
hole, i.e ~ 2 x 10 46 erg/s for a M sec = 1.4 x 10 8 solar mass 
secondary, equal to the total (beamed) luminosity of OJ287. 
This radiation does not arise from accretion of matter to 
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Figure 4. The difference energy spectrum of OJ287 between the 
outburst state and the pre-out burst state in 2005, with two-sigma 
errorbars. The lower points (crosses) have been corrected only for 
the Galactic extinction (hydrogen column density 3 X 10 20 cm - 2 ). 
The upper points (filled squares) include also a correction for the 
extinction in the OJ287 system, corresponding to the hydrogen 
column density of 9.3 X 10 20 cm~ 2 . Together they add up to the av- 
erage of the two measurements by Gosh & Soundararajaperumal 

(1995) . A Bremsstrahlung spectral curve for 3 X 10 5 K is shown as 
a dotted line, with the chi-squred probability of agreement with 
data P = 0.7. The fit with the power-law spectrum F v ~ v~ 0,8 is 
shown by a solid line and it gives the corresponding probability of 
P = 0.05 which is generally considered to imply a significant devi- 
ation between the theoretical line and the data. The best power- 
law fit to the uncorrected points follows F v ~ j/ -1,2 (\ 2 with a 
probability of agreement with data P = 0.01. The frequency unit 
is 10 14 Hz and the unit of vF v is 10 — 10 ergcm~ 2 s _1 . 

the secondary black hole but rather from a bubble of hot gas 
which is not accreted by it. In Section 5 we will return to the 
question of how much radiation may arise from the accretion 
process to the M sec — 1.4 x 10 s solar mass secondary. A more 
detailed calculation was carried out by Lehto & Valtonen 

(1996) . Scaling the accretion rate to 10 -3 MEdd where MEdd 
is the critical accretion rate of the primary black hole, and 
using the scaling laws of Stella and Rosner (1984) for the a g 
disk model of Sakimoto and Coroniti (1981), we find that 
the expected peak luminosity is 10 46 M ' 45 erg/s, i.e. roughly 
one half of the observed total luminosity. 

This extra luminosity is concentrated rather strongly 
in optical-UV. The peak bremsstrahlung flux at the optical 
V-band becomes S ~ 4a~ - 8 M 0A2 M sec mjy where M sec is 
the mass of the impacting (secondary) black hole in units of 
1.4 x 10 8 solar mass. For a typical a g value from observations 
(Dubus, Hamery and Lasota 2004, Qiao and Liu 2009) of 
a g ~ 0.3, S ~ 10Af sec mJy. This may be compared with the 
observed S ~ 6 mjy which agrees well with the calculated 
value when we take account of the likely extinction (Figure 
4), and use our reference value for the secondary black hole 
mass of 1.4 x 10 s solar mass. 

The disc at the impact point is heated initially to 
~ 5 x 10 6 K. The bubble torn off from the disc has the initial 
optical depth r ~ 60a~ °' 85 M°' 3 M S 144 , using the previous 
scaling for M and M sec . The expansion of the bubble by a 
factor of r 4/7 ~ Wet" 0A9 M 0A7 M° e f is required in order to 
make r = 1 at the time of the outburst. In adiabatic expan- 
sion the radiation dominated gas becomes cooler, and ac- 



quires the temperature of T ~ 5 x 10 5 a°- 49 M-°' 17 M-°- 82 K. 
Our fit to the energy spectrum in Figure 4 is fully consistent 
with this interpretation with a g ~ 0.3, M ~ 10 -3 Mem and 
M aec ~ 1.4 x 10 8 solar mass. 

Thus the mass value determined from the orbit solution 
M aec ~ 1.4 X 10 8 solar mass agrees well with the measured 
properties of the optical outburst. Our next step is to ask 
how much bigger the primary mass has to be in order that 
it can guarantee the stability of the primary accretion disc. 



3 STABILITY OF THE DISC 

The stability of accretion discs has been studied by many 
authors, in particular for the purposes of modelling inter- 
acting binary stars (e.g. Fragner & Nelson 2010). Larwood 
et al. (1996) studied different inclinations and binary mass 
ratios while Liu & Shapiro (2010) consider a wide range of 
mass ratios, but only coplanar orbits. The model for OJ 287 
has a binary orbit which is close to perpendicular to the disc 
plane and which has a large mass ratio. The stability of such 
systems has not been studied previously. 

We have modeled the accretion disc by a system of 6800 
particles which gives an adequate resolution in the study 
of the loss of disc particles due to the binary action. The 
particles are placed initially in a single plane approximately 
perpendicular to the binary orbit. The orbit integration is 
started around year 1850 and continued until year 2050 using 
the latest OJ287 binary model parameters which include 
the effects of the binary spin (Valtonen 2010a). The orbital 
elements of the particles are calculated at each time step. 
In the astrophysical disc model of Lehto & Valtonen (1996), 
the disc thickness is such that orbits with inclinations less 
than 6° remain in the disc (see Figure 4 of Valtonen 2007 
for the standard disc profile). 

We varied the secondary mass mi in order to get a 
set of runs with different mass ratios of binary black hole 
components, and for each run we calculate the fraction / 
of particles lost from the disc, i.e. particles with orbital in- 
clinations greater that 6°. Then the lifetime of the disc is 
determined from tuf e =200yr / f . The results are illustrated 
in Figure 5. 

The disc lifetime should be compared with the disc re- 
plenishment time which occurs in the viscous time scale 
(Ivanov et al. 1999): 

i 1ri 7 -0.8 iV— 0.4 1.4 /,% 

tvisc ~ 10 a M r yr (1) 

where a is the viscosity parameter, M is the accretion 
rate in units of O.OlMEdd, MEdd is the rate which gives the 
Eddington luminosity for the M = 1.84 x 10 10 solar mass 
central black hole, and r is the outer radius of the disc which 
has been perturbed, in units of 20 Schwarzschild radii. The 
parameters of OJ287 are such that £„ise ~ 10 7 yr. 

We notice that the disc particles are lost quite fast un- 
less the mass ratio M r is at least one hundred. The typical 
replenishment time scale of the disc by the viscous accre- 
tion process, tuf e ~ 10 7 yr, corresponds to the stability 
ratio ~ 130. This makes the minimum value of the primary 
mass ~ 1.8 x 10 10 solar mass which is quite close to the 
mass determined from the orbit, 1.84 x 10 10 solar mass. Liu 
& Shapiro (2010) find that the instability gap develops in 
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Figure 5. The stability time t;;^ e of an accretion disc when the 
mass ratio M r is varied. 



the disk when the mass ratio M r = 5 x 10 2 or less. The dif- 
ference is in the expected direction since tidal effects work 
less efficiently at high inclination than in the coplanar case 
(fvanov et al. 1999). 



4 BLACK HOLE MASS - BULGE 
LUMINOSITY CORRELATION 

Wright et al. (1998) measured the absolute magnitude of the 
bulge of the OJ287 host galaxy as M K = -28.9±0.6. In Fig- 
ure 6 we plot OJ287 in the black hole mass - K-magnitude 
diagramme using this value of Mk and the dynamical mass 
from Valtonen et al. (2010b). The other points are from a 
recent sample of Kormendy & Bender (2011a). The line of 
regression as well as the one sigma parallels are also drawn. 
We note that the point representing OJ287, with the dynam- 
ical mass from Valtonen et al. (2010b) and the K-magnitude 
from Wright et al. (1998) falls exactly on the regression line. 

Wright et al. (1998) derived also a lower limit for the 
host galaxy magnitude in the R-band and found that it 
is consistent with the R-band measurements by Hutchings 
(1987), Heidt et al. (1999) and Pursimo et al. (2002). Wright 
et al. (1998) point out that the K-band and R-band mag- 
nitudes imply about two magnitudes of extinction in the 
OJ287 host galaxy in the R-band, and larger hydrogen col- 
umn density that what we have discussed above. However, 
there does not need to be an exact correspondense between 
the measurements from the quasar and from the outer parts 
of the host galaxy since the lines of sight are different. Also 
there are large variations in extinction laws between differ- 
ent galaxies (Falco et al. 1999). It is quite likely that OJ287 
is a recent merger considering its double black hole nature 
(Iwasawa et al. 2010), and this impression is given also by 
the perturbed image of the host galaxy (Benitez et al. 1996) . 
Thus about two magnitudes of extinction in the R-band is 
not out of question since mergers tend to be dusty (Sanders 
et al. 1988). New observations of the OJ287 host galaxy in 
different colours would clarify the situation. 
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Figure 6. The black hole mass - K-magnitude correlation in the 
sample of Kormendy & Bender (2011), squares, and the OJ287 
point with magnitude error bars. The lines are the best fit re- 
gression to the Kormendy & Bender sample and the two sigma 
deviations from it. 



5 VARIABILITY STUDIES 

Studies of short time scale variability have been used to 
set limits on the black hole masses in OJ287 and in other 
quasars (Gupta et al. 2012). The idea is that the shortest 
periodic variability component may be related to the period 
of the innermost stable orbit of the radiating gas falling into 
the black hole. The same time scale, but shortened by the 
relativistic Doppler boosting factor, may also appear in the 
jet emission (Sagar et al. 2004). 

In case of OJ287 we have two time scales to consider: 
the scale related to the primary and the scale related to the 
secondary. Why should we consider the light output from 
the secondary whose mass is expected to be two orders of 
magnitude smaller than the mass of the primary? So far we 
have not resolved the two components separately but only 
record the total output from the two sources. Sundelius et 
al. (1997) calculated the infall rates of matter to both black 
holes, and found that the peak rates were about the same. 
The expected flux depends also on other factors such as 
the mass and the spin of the black hole. If the spin of the 
secondary is close to maximal, as would be natural after 
repeated accretion events from the constant high angular 
momentum source of gas (the primary disc), its luminosity 
may be an order of magnitude greater than the correspond- 
ing luminosity for a slowly rotating black hole. And since the 
primary black hole is likely to spin rather slowly (Valtonen et 
al. 2010a), the flux from the secondary could easily amount 
to about 10% of the flux of the primary. As to the Doppler 
boosting, it could be significant also in the secondary, since 
its spin axis is likely to be parallel to the angular momentum 
vector of its source of accreting gas, the primary disk. Thus 
the jet of the secondary could be beamed towards us just as 
the jet of the primary is (Valtonen and Wiik 2012). 

To find evidence for the emission of the secondary, we 
need to look at short time scale variability in OJ287. It 
has been found to be variable from 15 minutes time scale 
upwards on many occasions (Valtaoja et al. 1985, Dultzin- 
Hacyan et al. 1997, Gupta et al. 2012), and on one occasion 
the light curve has shown sinusoidal variations of the period 
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Figure 7. Observations of OJ287 brightness in the spring of 2005 
(squares), including the pre-outburst point in April 12 (crossed 
circle), compared with a 50 day periodic cyclic pattern. 



of 228 minutes (Sagar et al. 2004). If this variation is as- 
sociated with the last stable orbit of a maximally rotating 
black hole, the mass of the black hole is 1.46 x 10 s solar 
mass (Gupta et al. 2012), i.e. identical to the mass obtained 
from the orbit solution (Valtonen et al. 2010b). On the other 
hand, when the same frequency of variations is transmitted 
to the jet and is observed, the associated time scale is shorter 
by a factor of about 18.7, the Lorentz factor of the jet in 
the optical region. This is consistent with shortest observed 
variability time scale of 15.7 minutes in OJ287. 

Neronov and Vovk (2011) associate the 2009 October- 
November gamma- ray flare of O J287 also with the secondary 
black hole. 

As to the primary mass, the estimated period of the last 
stable orbit is around 100 days, the exact value depending 
on the spin of the black hole. The perturbation of matter 
from the last stable orbit, leading to its transfer to the jet, 
is affected by the tidal force of the companion. Since the 
tidal force is symmetric relative to the central black hole, 
the tidal transfer period from the last stable orbit should be 
one half of this value, i.e. about 50 days. Wu et al. (2006) 
reported the periodicity of 46 days in OJ287 during one 
observing season in the spring of 2005. The analysis of the 
14 year stretch of optical monitoring data from 1996 to 2010 
supports this finding (Valtonen et al. 2011). 

We have added 41 points from the 2005-2010 optical 
monitoring campaign of OJ287 (Valtonen et al. 2009) to the 
61 points of Wu et al. (2006) for the spring 2005; they are 
plotted in Figure 7. There is a clear variation at about 50 day 
period in this segment. A similar periodic variation extends 
up to January 2008; the peaks at 2005.85 and 2006.0 are 
the high points of this cyclic variation (Figure 1). The April 

2006 outburst and the three highest outburst peaks in late 

2007 and early 2008 (excluding the 2007 September impact 
outburst) are also part of the same 50 day cycle (see Figure 
3 of Valtonen et al. 2009). Within statistical uncertainties, it 
is possible that the period of one half of the last stable orbit 
has been detected, and that the primary mass is confirmed 
also in this way. 



6 DISCUSSION 

We find that the mass values for the two black holes in 
the OJ287 system, as derived from the orbital solution, are 
quite consistent with values obtained from astrophysical ar- 
guments. The chain of arguments starts from identifying the 
2005 outburst radiation as bremsstrahlung. Since its likely 
origin is the collision of the secondary on the accretion disc 
of the primary, this allows the determination of the mass of 
the secondary. The requirement that the disc is stable in- 
spite of the binary action puts a lower limit on the mass of 
the primary. It is very near the mass determined from the 
orbit solution. Thus we can have increased confidence in the 
correctness of the orbit solution. 

We place OJ287 in the black hole mass vs. K-magnitude 
diagramme and find that it lies almost exactly on the mean 
line determined from the lower black hole masses. This gives 
us confidence that this relation can be extended to up to the 
largest measured black hole masses. 

As a side result we find that the host galaxy of OJ287 
has a moderate amount of extinction which corresponds to 
the hydrogen column density of ~ 9.3 x 10 20 cm~ 2 , similar 
to values measured in some other galaxies from the 2175 A 
extinction feature (Wang et al. 2004). The dust and the per- 
turbed appearance of the host both suggest that OJ287 has 
acquired its second black hole in a merger of two galaxies, 
possibly leading to binary black hole evolution as calculated 
by Iwasawa et al. (2010). 

The observed periodicities in the light curves of OJ287 
are consistent with a model where both the slowly spinning 
primary black hole and the fast spinning secondary black 
hole contribute to the optical light variations (Neronov and 
Vovk 2011). In future it may become possible to see the 
orbital motion directly when the required microarcsecond 
resolution is achieved. It may also be possible to detect the 
motion of the secondary from its spectral lines. The orbital 
speed of the secondary is so fast that its spectral lines are 
shifted very far from their normal positions, outside the 
spectral region where the line searches have been carried 
out (Nilsson et al. 2010). 
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